Intensity Distribution of Modes in Surface Corrugated Waveguides 
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Exact calculations of transmission and reflection coefficients in surface randomly corrugated op- 
tical waveguides are presented. As the length of the corrugated part of the waveguide increases, 
there is a strong preference to forward coupling through the lowest mode. An oscillating behavior of 
the enhanced backscattering as a function of the wavelength is predicted. Although the transport is 
strongly non isotropic, the analysis of the probability distributions of the transmitted waves confirms 
in this configuration distributions predicted by Random Matrix Theory for volume disorder. 



PACS numbers: 42.25.Bs, 05.60.+W, 41.20.Jb, 84.40.Az 



Defects and inhomogeneities in waveguides are a sub- 
ject of increasing research for its influence on the propa- 
gation of both classical and quantum waves Q . Volume 
defects have recently been studied in connection with the 
propagation of electrons through two-dimensional (2D) 
wires |^ -H|. On the other hand, slight surface roughness 
has been considered both in electron transport through 
thin films as well as in optical fibers and waveguides 
P|, mainly in connection with attenuation. Recently, 
some analogies between electron transport and propa- 
gation of diffuse light through surface corrugated waveg- 
uides have been put forward H. However, the assump- 
tion of diffuse incoming waves limits the study to those 
systems in which the incident modes are mutually inco- 
herent. 

In this letter we address the more practical situation in 
light propagation where a single incident mode couples 
with different forward and backward modes As we 
shall see, the scattering from rough, perfectly conduct- 
ing, walls (see Fig.l) induces a strong and rich coupling 
to forward and backward modes. The coupling charac- 
teristics present significant differences with respect to the 
coupling induced by volume defects. For any incoming 
mode, and when the length of the system is larger than 
the localization length £, there is a strong preference for 
the forward propagation through the lowest mode. The 
coupling to backward modes presents a very interesting 
behavior, depending on the incident mode. The "exter- 
nal" modes (defined as those propagating modes with 
either the smallest or the largest transversal momentum) 
present enhanced backscattering factor larger than the 



others. In the case of the lowest mode, this factor ex- 
hibits remarkable oscillations as a function of the wave- 
length. We shall also present an extensive analysis of the 
statistical properties of the different transmission coeffi- 
cients. As far as we know, this has been carried out only 
for the case of volume defects. In the diffusive regime, 
the theoretical predictions of Ref. |9| (and subsequently, 
Ref. JTojl ) have recently been confirmed experimentally 
pd[ . A crossover from either the Rayleigh and Gaussian 
statistics in the diffusive regime to lognormal statistics 
in the localized regime, expected from numerical simula- 
tions jl2| , has been analytically deduced in Ref. jl3| . Al- 
though the mode coupling is quite different in the case of 
surface roughness, our calculated distributions of trans- 
mitivitties are fully consistent with the predicted behav- 
ior for bulk disordered waveguides. 

The corrugated part of the waveguide, of total length 
L and perfectly reflecting walls, is composed of n slices 
of length I. The width of each slice has random values 
uniformly distributed between Wq — 5 and Wq + S about a 
mean value Wq. We shall take Wq/8 = 7 and 1/5 = 3/2. 
The main transport properties do not depend on the par- 
ticular choice of these parameters, however. We consider 
s-polarized waves with the electric field parallel to the 
surface grooves (TE modes). Transmission and reflec- 
tion coefficients are exactly calculated by solving the 2D 
wave equation by mode matching at each slice, together 
with a generalized scattering matrix technique |7| |l^ , |l5| ] . 

For a given incoming mode i, the transmission (TV,-) 
and reflection (Rij) coefficients are defined by 
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qFotw b e j n g the total flux transmitted into the forward 
outgoing mode j, $^ acfc the total flux reflected into the 
backward outgoing mode j, and 3>|™ the total flux of the 
incident mode. The total transmission for the incoming 
mode i, is given by Tj = T^-, and G — £\ Tj is the 
normalized transmittance (conductance). Ensemble av- 
erages, denoted by (.), are performed over a thousand 
realizations of the corrugated waveguide (unless other- 
wise stated). 

The analysis of the conductance of surface corrugated 
waveguides is qualitatively similar to the electron trans- 
port in disordered nanowires. For short lengths the trans- 
port regime is diffusive with an effective mean free path 

l = (Nd(l/G)/dL)~ 1 , where N is the number of 
propagating modes. The transport regime changes as 
L becomes of the order of the localization length |H| 
£ = -(0(ln(G)>/0£) -1 , with £ = N£ within the nu- 
merical accuracy H. However, the analysis of the trans- 
mission coefficients reveals qualitative differences be- 
tween random surface and volume scattering. 

As a typical example, let us consider that the ra- 
tio between the waveguide width and the wavelength is 
Wo/A = 2.6. This allows five propagating modes. In Fig. 

1 we plot (ln(Tjj :)), as a function of the normalized length 
L/Wq for two different incident modes i = 1 and i = 3 
(Fig. 1 shows the average over 100 different realizations 
for each length L). All curves present the same linear 
dependence at large L. The inverse of the slope defines a 
length £ij = — (d(ln(Tij)}/dL)~ which does not depend 
on the incoming or outgoing modes (£r//Wo = 34.4 ±0.7 
for the structure calculated in Fig. 1). Exactly the same 
dependence is found for (lnTj). As a general result, and 
within the numerical accuracy, £y = = £ = Nt. 

In contrast with the expected isotropic mode distri- 
bution induced by scattering from volume defects, the 
calculation of and T$ shows a clear non isotropic dis- 
tribution, even in the diffusive regime (l < L < £) as 
shown in Fig. 2(a-c). In the localized regime (L > £), 
for any incoming mode i there is a strong preference to 
couple with the lowest transmitted mode j = 1. The 
coupling to forward modes decreases as their transversal 
(longitudinal) momenta increase (decrease). 

Concerning backward modes, Fig. 1 shows the aver- 
aged reflection coefficients (Rij ) as functions of the length 
L. Enhanced backscattering effects (EB) can be clearly 
seen in this figure. Coherent EB arising in the light scat- 
tering both from dense media p7| and rough surfaces 
Jl8[ has been intensively investigated For a waveg- 
uide with N propagating modes, Random Matrix Theory 
(RMT) predicts |l||2(J that the EB peak is 2 the scat- 
tering to any other channel, irrespective of the incoming 
mode. However, although not explicitly discussed, nu- 
merical calculations in waveguides with volume disorder 



B show factors other than 2 for EB of "external" modes. 
In order to quantify the EB peak, we define, for a given 
incoming mode i, an enhanced backscattering factor r){ as, 



m = (N-i) 



(2) 



Fig. 2d shows r\i versus L. While the "central" modes 
(i = 2,3,4) present a factor t) « 2, the "external" (i — 
1,5) modes have a much larger enhanced backscattering 
factor (ps 4). This qualitative behavior, which also con- 
firms indications for electron propagation in disordered 
wires [Q , depends on the wavelength. The backscattering 
factor for the lowest mode 771 oscillates above r\ w 2 as 
the wavelength A decreases (i.e. Wo/ A increases) having 
its maxima close to half integers of Wo/ A, which corre- 
spond to the appearance of a new propagating mode in 
the waveguide pi| . This behavior is highly correlated 
with the oscillations both in the mean free path I and 
in the localization length £ (tJ (the maxima in 77 corre- 
spond to minima in £). Just after the appearance of a 
new mode the backscattering dominates at this new ap- 
pearing channel. This gives a very large rj factor at the 
onset of the new mode, which slowly decreases to n = 2 
as A decreases. By contrast, the "central" modes present 
an almost constant enhanced backscattering factor near 
2. 

Let us now discuss the intensity probability distribu- 
tions of the transmitted waves. As shown above, the be- 
havior of the different averaged transmitivities for a sur- 
face disordered waveguide is qualitatively different from 
that obtained for volume disordered wires. Therefore, 
the question now is: does surface disorder change the 
probability distributions from those expected in bulk dis- 
ordered systems?. In the case of random wires, diagram- 
matic techniques combined with RMT pjl0|l show that in 
the diffusive regime the speckle pattern probability den- 
sity P(Tij) follows a Rayleigh statistics with stretched 
exponential tails, while the distribution P(Ti) is Gaus- 
sian with tail deviations. Nonperturbative calculations 
done in absence of time reversal symmetry predict that 
these distributions evolve into the same lognormal distri- 
bution as the length L increases beyond the localization 
length £ as expected from the numerical simulations of 
Ref. @. 

Fig. 3a contains the distribution P(Tj/(Tj)) for dif- 
ferent values of (Ti). Although for a fixed length, the 
distributions vary from mode to mode, we find that, as 
long as L < £, the distribution does not depend on the 
incoming mode i neither on the length L, but only de- 
pends on the averaged transmission coefficients. The in- 
set in Fig. 3a shows P(Ti/(Ti)) for i — 3,4, 5 all having 
the same average, (Ti) = 0.16, but corresponding to dif- 
ferent lengths of the disordered region (see Fig. 2c). For 
the same average, modes i = 1, 2 are localized and have 
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a different distribution (see below). The speckle distri- 
bution P(Tij/(Tij)), in the diffusive regime follows an 
almost perfect Rayleigh statistics P(x) = exp(— x) (see 
Fig. 3b). Clear deviations from this exponential distri- 
bution are observed for L < £ (not shown in the figure) 
and for L > £. An interesting point is that the depen- 
dence of (Tij) with length does not fix, by itself, the 
transport regime. While for I < L < £ all coefficients fol- 
low Rayleigh statistics, their dependence with L change 
completely from one to another as can be observed in 
Fig. 2a,b. 

Except for either small lengths L < £ or large trans- 
mitivities ((Ti) > 1/2), the qualitative behavior of both 
P{T t /(T t )) and P(T i:j /(T tj )) is basically the same as that 
obtained for disordered wires ||,[lo|,|l3| (compare our Fig. 
3 with Fig. 1 of reference ||l3|| ). Moreover, both distribu- 
tions evolve into the same lognormal distribution as the 
length L increases beyond the localization length. This 
is shown in Fig. 4 where P(ln(r)) is plotted versus hi(T), 
t = Ti,Tij, for different values of (ln(r)). Solid lines 
fitted there represent the lognormal distribution 
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with x = ln(T j; ),ln(T 4J ). and var(x) = a 2 (x) = (3/2)|(x)|. 
In the strong localization limit (L 3> £ ), where | (ln(r)) | ss 
L/f = L/(N£), we obtain a 2 {x) « (3/2)L/£. Neverthe- 
less, although the dependence is the same as expected 
from RMT Q, the variance of the calculated distribu- 
tion is 3/2 the mean in contrast with the expected value 
of twice the mean (cf. (11) of reference |13|). However, 
the analysis of the conductance G in the localized regime 
shows a lognormal distribution with a mean which is half 
the variance, in perfect agreement with RMT predictions 
[E^gol. This is shown in the inset of Fig. 4. The solid 
lines show the corresponding lognormal distribution (eq. 
3 with x = ln(G)) with var(x) = ct? ) = 2|(x)|. 

In conclusion, we have analyzed the coupling to for- 
ward and backward modes in surface corrugated waveg- 
uides. The coupling with backward modes yields en- 
hanced backscattering peaks, which can be larger than 
those theoretically predicted for the case of the lowest 
and highest incident modes. We have shown important 
differences in the transport coefficients with respect to 
the case of volume scattering. However, the probabil- 
ity distributions are the same as those obtained for bulk 
disordered waveguides and wires. In particular we have 
found a crossover from Rayleigh and Gaussian distribu- 
tions in the diffusive regime to the same lognormal dis- 
tribution in the localized regime in good agreement with 
the predictions of Random Matrix Theory. 
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FIG. 1. Averaged transmission (Tij) and reflection (Rij) 
coefficients (incoming mode i; outgoing mode j) as a func- 
tion of the normalized length of the disordered region L/Wo- 
Inset: Schematic view of the system under consideration. 



FIG. 2. (a) Averaged transmission coefficients (Tij) versus 
L/Wo (incoming mode i = 1, outgoing mode j). (b) The same 
as (a) for i=3. (c) Averaged transmission coefficients (Ti) ver- 
sus L/Wo. At (Ti) — 0.16 (horizontal dashed line) incoming 
modes i = 3, 4, 5 have the same probability densities as shown 
in the inset of Fig. 3. (d) Enhanced backscattering factor iji 
versus of L/Wo- Vertical long-dashed and dotted-dashed lines 
are drawn at L = £ and L — £ respectively. 

FIG. 3. (a) Distribution of x = Ti/ (Ti) for different values 
of (Ti). The inset shows the distribution for i = 3, 4, 5 having 
the same average (Ti) ~ 0.16 (see Fig. 2c). The smoothed 
curves are superimposed to the histograms to guide the eye. 
(b) Distribution of x = Tij / (Tij) for different lengths L/£. 
The straight line corresponds to the exponential distribution. 

FIG. 4. Distribution of ln(r) with r = T U T i3 for different 
values of (ln(r)). The inset shows the distribution of the con- 
ductance G for two different lengths of the disordered region. 
Continuous lines are fittings to lognormal distributions having 
a\ x) = 3/2(2;}, for x = Mr), and a\ x) = 2(x) for x = ln(G). 
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